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Abstract 



We study supersymmetric models without R parity and with universal soft 
supersymmetry breaking terms. We show that as a result of the renormaliza- 
tion group flow of the parameters, a misalignment between the directions in 
field space of the down-type Higgs vacuum expectation value and of the \x 
term is always generated. This misalignment induces a mixing between the 
neutrinos and the neutralinos, resulting in one massive neutrino. By means 
of a simple approximate analytical expression, we study the dependence on 
the different parameters that contribute to the misalignment and to m u . In 
large part of the parameter space this effect dominates over the standard 
one-loop contributions to m u ; we estimate 1 MeV < m u < 1 GeV. Labo- 
ratory, cosmological and astrophysical constraints imply m v < 100 eV. To be 
phenomenologically viable, these models must be supplemented with some 
additional mechanism to ensure approximate alignment and to suppress m u . 
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I. INTRODUCTION 



The field content of the Standard Model (SM) together with the requirement of SU (2) L x 
U(1)y gauge invariance, implies that at the renormalizable level the most general Lagrangian 
possesses additional accidental £7(1) symmetries. The U(l) generators correspond to Baryon 
(B) and Lepton flavor (Li) charges. The conservation of B , Li and hence of total Lepton 
number (L = J2iLi) naturally explains nucleon stability as well as the non observation of 
L and Li violating transitions. This nice feature of the SM is lost in its Supersymmetric 
(SUSY) extensions: once the SM fields are promoted to superfields, additional gauge and 
Lorentz invariant terms are allowed, which violate B, Li and L at the renormalizable level. 
To forbid these dangerous terms, a parity quantum number R = (—\ s j 3B + L + 2S (S being 
the spin) is assigned to each component field, and invariance under R transformation is 
imposed. Therefore, in SUSY frameworks B and L quantum numbers are assigned ad hoc 
to the superfields to reproduce the accidental symmetries of the SM, and ensure the absence 
of non observed transitions. However, R parity is by no means the only symmetry which 
allows for building viable SUSY extensions of the SM |L|. From a phenomenological point 
of view, the first priority is to ensure the absence of operators leading to fast nucleon decay, 
and in this respect other discrete symmetries can be more effective than R. This is because 
R parity only forbids dimension 4 B and L violating terms, but does not forbid dimension 
5 operators which can still be dangerous, even when suppressed by factors as large as the 
Planck mass. Some interesting alternatives to R parity exist, which forbid dimension 4 and 
5 B violating terms but do not imply the same for the L non-conserving terms QXJ , and thus 
imply a rather different phenomenology from models with i?-parity. Since a mild violation of 
L can be phenomenologically tolerated, SUSY extensions of the SM with highly suppressed 
B violation but without R parity and without L number, represent interesting alternatives 
to the Minimal Supersymmetric Standard Model (MSSM). Two new types of Lagrangian 
terms characterize this class of models: (i) renormalizable interactions responsible for L and 
Li violating transitions; (ii) superrenormalizable terms which mix the three lepton doublets 
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with the down-type Higgs. These mixing terms are present because the four hypercharge 
Y = — 1 doublets transform in the same way under the full gauge group. They also imply that 
after the gauge symmetry is spontaneously broken, all the color singlet fermion fields with the 
same electric charge are mixed (left and right handed charged leptons mix with Higgsinos 
and winos, neutrinos mix with neutral Higgsinos and with the zino). This situation can 
appear phenomenologically untenable, however the mixing acquires a well defined physical 
meaning only when a physical basis for the various fields is defined. We define the down- 
type Higgs as the particular combination of the four Y = —1 doublets which acquires a non 
vanishing vacuum expectation value. If all the superrenormalizable terms in the Lagrangian 
are such that in this basis the remaining three combinations are decoupled from the Higgs, 
we can still assign to the fields a lepton number which is violated only by the renormalizable 
interactions (i) while, at lowest order, it is conserved by the mass terms and by the gauge 
interactions. In the class of models where the soft terms responsible for SUSY breaking 
are universal, the conditions required to realize this scenario seem to be satisfied. Since 
minimal SUSY extensions of the SM generally belong to this class, most of the literature 
on SUSY without R parity concentrated in studying the effects of the interaction terms 
(i) [^-|3T|, while less attention has been payed to the consequences of (ii) . However, 

even in the minimal models, universality of the soft terms holds only at some high energy 
scale where these terms are generated. The set of loop corrections induced by the terms 
(i) imply deviations from universality for the low energy parameters, and this unavoidably 
results in the appearance of the terms (ii) which therefore have to be always included in the 
i?-parity nonconserving superpotential. Moreover, since the scale where universality holds 
can be as large as the Planck scale, deviation from universality at low energy can be relevant 
and imply that the effects of the Renormalization Group (RG) induced superrenormalizable 
terms cannot be neglected. 

In Section II we first present a qualitative discussion, based on symmetry considerations, 
of the mass spectrum for the color singlet fermions. We also review the conditions for the 
alignment in field space between the down-type Higgs vacuum expectation value and the 
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[i term |§[T4| which plays a crucial role in ensuring the suppression of neutrino masses. In 



Section III we derive a simple formula which parametrizes the RG induced misalignment, and 
we discuss the main dependence of this effect on the model parameters. The fermion mass 
spectrum is discussed quantitatively in Section IV, where the mass of the heaviest neutrino 
arising from the misalignment is estimated. Section V contains a brief review of the main 
laboratory, cosmological and astrophysical constraints on the neutrino mass, which can be 
translated into constraints on the relevant parameters responsible for the misalignment. Our 
results are summarized in Section VI. 

In many aspects our analysis complements recent works that discuss the same effect 
32H34| . Refs. [ 32] and |34]] restrict their analysis to models in which L violation enters only 



through the bilinear terms (ii). The renormalizable interactions (i) arise from the Yukawa 
terms only after the fields are rotated to a basis where the lepton doublets are decoupled 
from the Higgs. Therefore the form of these terms is not general, but is determined by 
the corresponding Yukawa couplings with a proportionality factor accounting for the field 



rotation. Issues analogous to the ones studied here are also addressed in Section IV of |B3 
where some results corresponding to specific choices of the parameters are presented. A brief 
discussion of these effects is also given in Section VI of [|T3|]. Most of the results presented in 
these studies are given in numerical or graphical forms, which render difficult to appreciate 
the details of the physics involved. In the absence of analytical results it also appears 
awkward the task of taking properly into account these effects in future studies of SUSY 
models without i?-parity. In the present work we study the general .R-parity violating case by 
including all the terms consistent with the SM gauge symmetry and with B conservation. 
In contrast to previous works, our approach is essentially analytical. We give a simple 
basis independent expression for the RG induced misalignment which highlights its physical 
meaning. We present an analytical formula for the neutrino mass that shows explicitly the 
main effects involved, and makes it easy to appreciate the various interrelations between 
the different parameters of the model. All our main results are summarized in a few simple 
expressions that can be easily used for investigating further this class of models. 



II. ALIGNMENT 



In this section we examine the qualitative features of the fermion mass spectrum which 
can be expected in SUSY models without R parity. The SUSY extension of the SM contains 
eight color-singlets chiral multiplets, corresponding to the up-type Higgs field, three right- 
handed leptons, three left-handed leptons and the down-type Higgs doublets. Under the 
electroweak gauge group SU{2) L x U(l)y their quantum number assignments are 

H u ~ (2, 1), 

I ~ (1,2), (i = l,2,3) 

H a ~ (2, -1), (« = 0,1,2,3). (1) 

Here H a denotes collectively the supermultiplets containing the down-type Higgs and left- 
handed lepton doublets, which in the MSSM are distinguished by different R-parity assign- 
ments. If R-parity is not imposed, the gauge interactions posses a global SU(4) symmetry 
corresponding to rotations of the four H a superfields ||34|| . However, other terms are gener- 
ally present which select some preferred directions in SU(4) field space. The relevant terms 
that break the symmetry in the fermion sector are: 

(a) The bilinear superpotential term 

/i a H a H u , (2) 

provides a mass for the fermionic component of one combination of the H a doublets 
(the Higgsino). The symmetry is broken down to SU(3) acting on the three combina- 
tions orthogonal to the Higgsino. 
(6) The vacuum expectation values (vevs) 

(H a ) = v a , (3) 

which contribute to the spontaneous breaking of the electroweak (EW) symmetry, 
induce a mixing between the neutral members of the H a fermion doublets and the 
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neutral gauginos, thus breaking the symmetry down to SU(2). A second combination 
of the neutral members in H a acquires a mass in this way. Since the vector v a fixes 
a direction only for the neutral fields, no additional charged fermion becomes massive 
at this stage. 

(c) Finally, the following trilinear terms in the superpotential break the symmetry com- 
pletely: 

X a /3k H a Hp £ k + \' a - k H a Qj dk ■ (4) 

Here Qj and dk denote the quark doublet and down-quark singlet superfields, and 
X a /3k = ~^f3ak due to the antisymmetry in the SU (2) indices. For the charged fields 
the breaking is induced at tree level by the X a pk couplings, which generate three new 
vectors in SU(4) space corresponding to the mass terms {m) a k = \ a pk v Pi (k = 
1, 2, 3). In the neutral sector the residual SU(2) symmetry is broken only at the loop 
level, through quark-squark and lepton-slepton loop diagrams which generate the mass 



terms 15,251 
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Here m d is the (i-quark mass matrix which arises at tree level from the second term in 
(£§), Mui i s t ne left-right sector in the d-squark mass-squared matrix, M[ R is the left- 
right sector in the mass-squared matrix for the charged £j-H a scalars, m represents a 
slepton or squark mass, and the expression holds at leading order in M\ R j fh 2 . 
The qualitative features of the fermion mass spectrum for the fields in (P arising from 
the pattern (a)-(c) are the following: 

(1) Only one combination of the charged H a acquires a large mass of order fj, = (fiafia) 1 ^ 2 
(or of the order of the EW breaking scale) while the remaining three charged fermions 
get masses proportional to (arbitrarily small) Yukawa couplings. 

(2) Two neutral combinations of the H a doublets acquire small masses only at the loop 
level (c), while other two get large tree level masses as a consequence of (a) and (6). 



Since three neutral fermions (neutrinos) are known to be very light, (2) represents a major 
challenge for reconciling this scenario with our experimental knowledge of the fermion mass 
spectrum. The qualitative features of the predicted spectrum can be reconciled with the 
observations if some mechanism ensures that 

V a OC fl a . (6) 

If this relation is satisfied, no new direction is singled out by the vevs v a , SU(3) still 
remains a good symmetry after (b) and three neutral fermions acquire their mass only 
through the loop effects (c). As we will see below, in general (^) cannot be enforced as an 
exact (low-energy) relation. However, since the tree level mass which is induced at stage (6) 
is proportional to the amount of SU(3) breaking, an approximate alignment between v a and 
fi a can be sufficient to avoid conflicts with the limits on neutrino masses. 

The conditions for v a and \i a alignment were studied in |]3|JT^|. The direction of v a is 
determined by the minimum equations for the scalar potential, which depend on the soft 
SUSY breaking terms 

HaHp , (7) 

and on fj, a . Terms proportional to X a/ 3j and X' ai j as well as the soft SUSY breaking trilinear 
A terms which also carry SU(4) indices, always involve a charged field and hence at lowest 
order do not contribute to determine v a . Relation (^j) holds if the following two conditions 
are satisfied MM: 



(A) y. a is an eigenvector of m 2 HaHg : fh% aHp ftp = m 2 y. a ; 

(B) B a is proportional to \i a : B a = B \i a . 

To show this, let us rotate the H a fields to the basis (H\\,H±) where H\\ = fi a H a / \i and 
H± denotes the three combinations orthogonal to H\\ . According to (A), in this basis 
rfijj Hg = rh 2 5\\p {(3 = ||,_L) while (B) implies that, like fi a , also B a has the only non- 
vanishing component along H\\ . Then the solution of the minimum equations corresponds 
to (H±) = . The vector v a is thus aligned with jx a : Vd = (vaVa) 1 ^ 2 = (Hn) = \i a v a j 'fj, . 



III. RENORMALIZATION GROUP INDUCED MISALIGNMENT 



In models where SUSY breaking is induced by universal soft breaking terms, both condi- 
tions (A) and (B) of the previous section hold. However, these conditions are exactly satisfied 
only at the scale Au where the universal terms are induced. After defining B a = B a/3 Hp , 
the universality conditions at A^ read 



(B a /s)A u = B v 5 a p . (8) 

As a result of the RG running of the parameters, these relations become only approximate 
at low energy, and a misalignment between \x a and v a is generated. The deviations from 
conditions (A) and (B) at a generic energy scale can be parametrized as 

(g B al3 - S <*p = A f/3 Hp (9) 

where fh 2 ~ (det ™ 2 HaHl3 ) 1 ^ and B = (BaBa) 1 / 2 . To estimate the misalignment induced by 
the RG running, we write the result of the minimization of the (low energy) scalar potential 
as 

Va = -4 [ S a p + A Q/3 ] (ip , (10) 

A 4 

where A a p = A™^ + A^g accounts for the misalignment induced by violations of conditions 
(A) and (B), and the normalization factor fx ~ /i. 

It is now convenient to introduce two unit vectors e v a and with components v a /vd and 
Ha/ A* • The misalignment can be quantified by means of an angle £ defined as 

sin£= |e" Ae v \ . (11) 

By means of fllQl) we obtain 

sm 2 £=^E [K^-e^ M )e^] 2 

= e M ■ A 2 ■ e M - (e M ■ A • e M ) 2 . (12) 



We note that all factors proportional to 5 a p in A a! 3 cancel in fll2|). Therefore, in com- 
puting sin£ it is sufficient to retain only the terms which carry non trivial SU(4) indices. 
In particular, the contributions proportional to the up-quarks Yukawa couplings and to the 
gauge couplings can be dropped off, and only terms involving the couplings in (Q) and the 
corresponding trilinear soft SUSY breaking A terms need to be kept. An approximate an- 
alytical expression for A a/ 3 , obtained by assuming constant coefficients and by integrating 
the RG equations in one step, is derived in the Appendix. Proceeding in this way, it is 
possible to single out the main effects which generate misalignment and to keep track of the 
various interrelations among different parameters. Besides universality, in the following we 
will also assume that at Ay the trilinear soft SUSY breaking A terms are proportional to the 
corresponding couplings in This assumption is made only for reasons of simplicity, since 
it allows factoring out the overall scale Au of the soft-breaking terms and this gives simpler 
expressions. However, the results of the analysis do not depend on this assumption, and it 
is straightforward to replace terms like AyX^ with the more general parameters A a ^ . At 
the EW scale ~ mz , the RG induced misalignment matrix reads 

= ^3 ^3 + + (X al iXpji + SX'aijX'pij) , (13) 

where ty = log Mz / Ay , Ay , By and m v are the soft SUSY breaking parameters at Ajj , 
and only terms inducing 577(4) rotations are displayed. In (|13D the term proportional to 
AujBu comes from the running of B a p, while all the others originate from fri^g- We learn 
the following : 

(a) The RG induced v a -fi a misalignment originates mainly from the running of the the 
soft-breaking scalar masses. For Ajj ~ Bjj ~ rhjj this effect dominates by a factor 
of ~ 4. Only if By <C Ay <C my the misalignment is dominantly induced by the 
evolution of the B terms. 

(6) Apart from fine tuned cancelations and as long as L is a broken symmetry, there is no 
limit for the soft-breaking terms in which alignment can be recovered. 
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(c) If Ay <C By, mjj, the misalignment is independent of the initial values of the soft- 
breaking parameters. 

(d) Since SU(4) rotations in the evolution of B a p are induced only by terms proportional 
to Ay (see the Appendix) if Bjj = the third term in the first parenthesis in ( p~3|) 
is unity If Ajj = Bjj = the B term does not contribute (at this order) to the 
misalignment. 

From (|12D we see that sin£ is a basis independent physical parameter. It is convenient 
to give its explicit expression in a specific basis. We define the down-type Higgs field H d 
at the EW scale by the condition (H d ) = v d (that is H d = e v a H a ) and we choose the 
basis {H d ,Li} where Lj are three states orthogonal to H d . In this basis e v a = <5 0a , while 
hfj = X' a j k e v a and hi k = \ a pk^ v a (with h l Qk = h e ak e v a = for the antisymmetry of the A 
couplings) are the H d Yukawa couplings to the fermions. After inserting ( |13D in (|1^) and 
using e£ ~ (5 af3 - A a/3 )e^ ~ 5 0a we obtain 



sin 2 f ~ A m A i0 



tu \ 2 Afj , Ajj ^ 



3 + ^T + 

u D u 



^7r 2 / \ mh B. 



where it is understood that Xijk and A^- fc are now the couplings (§) rotated to the {H d ,Li} 
basis.Q From (]14 ) we learn the following: 

(e) To generate misalignment is enough to have at least one of the L violating A or A' 
couplings (or one of the corresponding A terms, if the assumption of proportionality 
is dropped) non-vanishing. 

(f) Assuming no particular suppression of the R-parity violating 6-quark couplings A^ 33 
with respect to the couplings involving the first two families, the dominant contribution 
to the misalignment is proportional to the 6-quark Yukawa coupling hf 3 . 



*For simplicity we do not distinguish between the couplings at the EW scale and at An- For our 
approximate solutions of the RG equations the difference is formally of higher order. 
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(g) Since only the leptons and <i-quarks couplings appear in flT4]), the misalignment de- 
pends strongly also on the value of tan/3 = v u jv& (where v u = (H u )). For the leading 
contributions we obtain 

2 / \ 



sin 2 £ 



3tu m b 



in 2 v 



A 



A 



(15) 



where v = (v% + f 2 ) 1 / 2 ~ 246 GeV. If in addition we assume that the L violating 
couplings are proportional to the corresponding Yukawa couplings (as is the case in 
models based on horizontal symmetries p^ , [15| ) then sin 2 £ ~ tan 4 /?. 



IV. FERMION MASS SPECTRUM 

In this section we investigate the consequences of the RG induced misalignment on the 
fermion mass spectrum. Numerically, the factor in square brackets in (|T5| ) is at most of 
order 1% resulting in sin£ <C 1 and approximate low energy alignment. As we will see, this 
implies that L violation in the mass terms is small, and that a distinction between 'leptons' 
and charginos and neutralinos is still a meaningful one. We define the 'right handed leptons' 
as the mass eigenstates having as main components the three SU(2) singlets li. Their mass 
partners are the 'left handed charged leptons' which are dominantly combinations of just 
the Y = — 1 doublets. Their neutral SU(2) partners constitute the main components of the 
neutrinos, while the remaining mass eigenstates are the charginos and the neutralinos. 

The mass matrix for the charged fermions M c is 5 x 5, with rows corresponding to 
{W~, H~}, and columns to {W + , H+, 1^} : 

/ M 2 -j=v u lx3 

Here 0i X 3 denotes a zero 1x3 block and M 2 is the SU(2)l gauginos Majorana mass. In the 
{W~ , Ha, L^} basis (in which we denote particles and superparticles according to the usual 
convention) this becomes 
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K 



\ 



(17) 



M 2 V2m w sin (5 0i x3 \ 

/icos£ 0i X 3 
3 xi /i(e"Ae")i /if^y 

where mw — gv/2. As expected the charged lepton masses originate from the Yukawa cou- 
plings to H d , and their mixings with the charginos, induced by the (3 2) block, is suppressed 
at least as sin £ . (If sin £ is not too small, this mixing could still give rise to interesting 
processes like Z -> W + L~ , LfLj , £f£j (i ^ j), etc.) 

The full neutralino mass matrix is 7x 7. In the basis with rows and columns corresponding 
to {B, W 3 , H°, Hi} it reads 



V 



Mi 





gta,n8 w 

2 ^ 


g tan 
2 







M 2 


_£ v 

2 ^« 


2 ^ a 




2 ^ 


_a v 

2 









gtan0 w 

2 Uq 


2 U ° 




04x4 


J 



(18) 



Here Mi is the U(l)y gaugino mass and 9w is the weak mixing angle. In the {Hd, L^} 
basis v a = (vd, 0,0,0) and \i a ~ (/i, —A i0 /i). M n gives 5 massive states and two massless 
ones. Four massive states correspond to the neutralinos while the fifth one, a neutrino, 
corresponds to a combination of the neutral members in Lj. The mass of the neutrino is 
given by 

det'M„ 



m v 



(19) 



(=0 



det'M n 

where det' denotes the product of the nonvanishing eigenvalues of the respective mass ma- 
trices. We have 



det'M n = /i 2 M^ m| cos 2 (3 sin 2 £ , 



det'M n 



5=0 



[x M^ m| sin 2(3 - /i 2 Mi M 2 , 



(20) 



where rriz = gv/(2cos9 w ) and My = Mi cos 2 6w + M 2 sm 2 w is the photino mass term. 
The final expression for the mass of the neutrino reads 



sin 2f3 



fxM 1 M 2 ~ 


-l 


'3tum b ' 


'( 


M^m 2 z 




87T 2 V 





At 



A 



i33 



(21) 
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From this expression we see that the numerical value of m u depends on several parameters, 
and this explains why it is not easy to derive any simple scaling behavior from a numerical 
study of these effects. However, the leading behaviors are as follows: if fxMilvr? z 3> sin 2/3 
(which is more easily satisfied in the tan/5 3> 1 limit) then m u oc m z /M2, practically 
independently of fx and tan/3. Therefore, as was noted in [32|, m u vanishes in the limit of 



very large SUSY breaking scales. In the opposite limit, which is consistent only for moderate 
values of tan/5 ( < 5) and for small values of Mi and fx , m„ is approximatively proportional 
to fx tan (3. It is interesting to note that in the first equation (|20| ) cos 2 (3 cancels against 
the explicit 1 + tan 2 /3 term in the misalignment parameter sin 2 £ fll5|), leaving only a mild 
dependence on tan/3 in the final result. However, if the .R-parity violating couplings are 
proportional to the Yukawa couplings, an implicit tan 2 (3 dependence from the misalignment 
is still present in the last parenthesis in fl2ip . We also note that the first square bracket in 
cannot approach zero, since is bounded by the lower limits on the neutralino masses. 
Banning possible fine tunings, for natural values of the parameters we obtain that fx divided 
by the first bracket yields a dimensionful factor ~ 10 — 100 GeV|] The square of the second 
bracket provides a suppressing factor in the range 10 _3 -10 -5 corresponding respectively to 
A {7 ~ m Planck and ~ 10 5 GeV, where the second value is typical of gauge mediated SUSY 
breaking scenarios |36|-^l|. Finally, the square of the term containing the soft breaking 
parameters yields approximatively a one order of magnitude enhancing factor. 

As a result, in the absence of further suppression from the A' couplings, we would estimate 
1 MeV < m v < 1 GeV. 

Before concluding this section, it is interesting to compare the RG induced effects 
on m u with the standard contributions to the neutrino mass matrix from one-loop dia- 
grams |p!5| , p5| -|3T|] . From ([5]) the corresponding leading term reads 



^The lower limit can be pushed down to 1 GeV in a phenomenologically viable scenario in which 
fx is very small resulting in two neutralinos not much heavier than a few GeV [35| 
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3 m 2 

<^^^( A -^^KA- (22) 

The misalignment yields the dominant contribution to the mass of the heaviest neutrino as 
long as 

,2 1 1/2 



tjj > 



57T . 

— F(m soft ) 



(23) 



where the dimensionless function F depends in a complicated way on the various soft SUSY 
breaking parameters, as well as on fi and tan/?. For most values of F fl23|) is satisfied as 
long as Ay > 10 5 GeV, and thus in general the induced misalignment gives the leading effect. 
This implies that predictions for the neutrino masses in models without i?-parity and with 



high energy alignment based solely on an estimate of the loop contributions |5T|, should be 
modified to include this effect. 

On the other hand, F is maximal when Ay <C %, By and in this limit we obtain 
F < (1/g 2 ) (fiM 2 /rh 2 ) tan f3. This situation is interesting since Ay = can arise in gauge 



mediated SUSY breaking models |j37H - In this case, for values of the relevant parameters 
such that F ~ (l/g 2 ) and for Ay < 10 6 GeV the two effects yield contributions which can be 
comparable in magnitude. Finally, for small Ay and tan (3 rather large (11M2 tan/3/m 2 > 25 ) 
the one-loop contributions (p2|) dominate over the misalignment effects up to Ay = mpi anc k , 
and thus determine the mass of the heaviest neutrino. 



V. PHENOMENOLOGICAL CONSTRAINTS 

As we have seen, in SUSY models without R-parity and without Lepton number the in- 
duced v a —fi a misalignment results in one massive neutrino which, in the absence of suppres- 
sion of the L violating trilinear couplings, is naturally in the range 1 MeV < m v < 1 GeV. 
In this section we argue that laboratory and cosmological constraints imply that this win- 
dow is excluded. In addition the massive neutrino of these models is very likely stable on a 
cosmological time scale, and thus the cosmological limit m u < 100 eV applies. As a conse- 
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quence, to render these models phenomenologically viable some mechanism to suppress the 
L- violating couplings in ( f2T|) down to J2i -^33-^33 ~ 10~ 4 - 10~ 7 is called for.f] 

The flavor composition of the neutrino is determined by the relative rotation in the 
{Ha, Li} basis arising from the diagonalization of the submatrix in (|U|) containing A,o (for 
the z/'s) and of the Yukawa couplings matrix h e h e (for the left-handed leptons). This can 
be studied only by specifying further the model. We will avoid doing this, and we will 
conservatively assume that our massive state is mainly z/ T , so that the laboratory limit on 
the mass is m u < 23MeV ||45|| . With regards to the neutrino mixing angles, the discussion 
below is purely phenomenological and does not require any theoretical estimate. 

Cosmological considerations of the age and the present energy density of the Universe 
provide constraints relating the mass and lifetime of the neutrino. For masses in the range 



100 eV - a few MeV the constraint reads m v T v < 2 x 10 MeV sec j46[ . When charged 
particles are present in the final state, a stronger bound from the absence of distortions in 
the cosmic microwave background radiation (CMBR) applies, r v < 10 4 sec. Detailed studies 
of the effects of a massive v T during the nucleosynthesis era imply an even stronger limit 
t v < 10 2 sec, for masses in the range 0.5 MeV - 35 MeV and independently of the decay 
modes P7^p^[. For visible decay modes (final states containing 7 or e ± ) also a lower bound 
on r v exists, r v > 10 8 sec. This bound follows from the limits on the gamma-ray fluence 



around the time when the neutrinos from the Supernova 1987A were detected ||53|| . This set 
of constraints already suggests that m v > 1 MeV is very likely ruled out. 



■^Such a suppression can be easily accommodated in models for fermion masses based on Abelian 
horizontal symmetries [^-44|. The required suppression for the case when the soft SUSY breaking 



terms are not universal and the misalignment arises as a tree level effect was studied in [14,15]. 
In contrast to that case which required horizontal charges for the Lj doublets larger than ~ 7 , in 
the present scenario a sufficiently small neutrino mass is obtained with the more natural values 
Q H {U) ~ - 3. 
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In order to avoid some (or all) of these constraints, the massive v should decay fast 
enough, and preferably into invisible final states. However, most likely the dominant v 
decay mode is v — > e + e~v e which proceeds via H^-mediated tree level diagrams. All other 
decay modes, as v — > 3z/£ or v — > ji/£ (£ = e, fi) are flavor changing neutral current (FCNC) 
processes, and are more suppressed. In particular, the invisible decay mode into three 
neutrinos involves the FCNC Zvvg vertex, which is quadratic in the neutrino mixing with 
the isotriplet neutralino W3, and hence very small. For the leading decay mode the lifetime 



is 



m, 



where ~ 2.2 x 10~ 6 has been used. On the other hand, peak and kink searches in 7r , K and 
(3 decays yield stringent upper limits on |t/i^| 2 . We have |54| \U eu \ 2 < 5 x 1CT 6 (1 x 10 -4 ) for 
a mass of about 20 MeV (5 MeV), implying lifetimes in conflict with the nucleosynthesis (and 
CMBR) constraint. For smaller masses the laboratory limits on the mixing parameters are 
less stringent. However, below about 3.5 MeV (1.5 MeV) the constraint from nucleosynthesis 
(and CMBR) is not satisfied even for maximal mixing. For m u < 1 MeV only FCNC decay 
channels are open, implying that also the weaker mass-lifetime constraint from the age of 
the Universe is not evaded. Therefore, we conclude that independently of the mass and 
mixing angles the v decay rate is not fast enough to evade all the constraints, and the limit 
for cosmologically stable neutrinos holds. 



VI. CONCLUSIONS 

In this paper we have presented an analysis of SUSY models without R parity and 
without Lepton number. We have shown that even when universality of the soft SUSY 
breaking terms is assumed, at low energy the vector v a of the vevs of the hypercharge — 1 
doublets H a is not aligned with the vector fi a of the generalized /z-term fi a H a H u . The 
misalignment is induced by the renormalization group flow of the parameters from the 
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scale where the soft SUSY breaking terms are generated, down to low energy. We have 
derived a simple analytical expression which describes the dependence of the misalignment 
on the relevant parameters. Our treatment is basis independent, and shows that this effect 
cannot be rotated away or neglected. In the basis where the fields are physical, the bilinear 
superpotential terms \XiLiH u which violate Lepton number by one unit, are always present. 
A major consequence of v a -/i a misalignment is that one neutrino becomes massive, and the 
mass induced in this way is generally larger than the contributions from one-loop diagrams. 
We have estimated that in the absence of additional suppression 1 MeV < m v < 1 GeV. 
A brief analysis of various laboratory, cosmological and astrophysical constraints strongly 
suggests that this neutrino is cosmologically stable, and thus its mass must be below 100 eV. 
This bound can be translated into a constraint on the .R-parity violating trilinear couplings 
J2i ^33X33 ~ 10 4 ~~ 10~ 7 . We conclude that, to be phenomenologically viable, SUSY models 
without -R-parity must be supplemented with some mechanism (as for example a horizontal 
flavor symmetry) yielding a sufficient suppression of these couplings. 
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APPENDIX A: 

In this Appendix we compute the misalignment matrix 

= A3, + A™ (Al) 

induced by the RG evolution of the soft SUSY breaking parameters rh 2 HaHf3 and B a/3 from 
the high scale down to the EW scale. The running is controlled by the RG equations 

drh H aHf3 m 2 - a 
dt 16tt 2 a/3 ' 
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dt ~ 16tt^' (A2) 
where we have factored out the overall scale fh 2 and A of the soft SUSY breaking masses and 
of the trilinear soft breaking terms. The boundary conditions at Ajj , where by assumption 
universality holds, are given in @. We solve (|A2|) in first approximation, neglecting the 
scale dependence of the coefficients. This gives 



^a(3 = ~FT~ B af3 ~ S a fi ~ — (G a/3 ) Au , (A3) 

where t v = log (M z /Au) . The SUSY RG equations including R-parity violation have been 



recently presented in a number of papers |33| , |55| , |56| . The equations for the soft SUSY 
breaking terms can be read off from |33]. For the running of rh 2 H H in (|A2|) we have 



m 2 G% = C H S < pHa + m|,- + C Q J m% Qj + C d ^ m% dj + A 2 C% + C G 5 a , (A4) 
where 

9 ^aai^/3pi + 3(\' aij \' pij 5f3 a + \' f3ij \' pij 5 aa ) 



syH pa 

U a/3 


X aryi " p-fi ®(3a 




fil ij 
U a(3 


2 \or/i^fryj 




s-yQ ij 

U a/3 


= 6 ^aik^fijk 




n d ij 
° a/3 


= 6 ^aki^(3kj 




°a/3 


2 (A Q ,'yjA^j 


+ 3 ^aij^(3ij) 



C G = - £ 9lY a mla ~ 19\M\ - 6g 2 2 Mi . (A5) 

cr=all 

The running of the B term is determined by 

G% = D O + D8 a p (A6) 

with 

D a p = 2 \ \ Q . 1 i\p 1 i + S^'aij^'pij) 

D = 2 (3t%h% + gf^ + 3^) , (A7) 



where are the up-quark Yukawa couplings. Using the boundary conditions @ together 
™ th = (™ 2 QiQj )Av = (™hMu = S H we obtain| 



:/3 Mi/ 



A 2 



6 + 2 -rj- (XoeyiXp^i + 3 A' Qifc A' 



rn 



/3ik) 



{Ga8)Av — 2 (AcryiA/^j + 3X' ai ■) + 2 



=all 



2 Mi/ 



mf 



J a/3 



3^ + ^ + 3^ 



2 M a 



(A8) 



A v "Act. 

where Mjj is the universal gaugino mass. The terms in square brackets which are pro- 
portional to 5 a/ 3 do not generate misalignment, and for our purposes can be dropped off. 
Inserting the relevant terms of equations ( |A8| ) into ( |A3| ) gives the expression ([T3|) for the 
misalignment matrix A a p. 



Hn models where SUSY breaking is communicated to the visible sector via gauge interactions, the 
SUSY breaking terms are flavor-symmetric but not universal. The corresponding modifications of 
equations ( |A8[) are straightforward, and do not affect our conclusions. 
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